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Abstract 

In this paper Janus black holes in AdSs are considered. These are static solutions of an 
Einstein-scalar system with broken translation symmetry along the horizon. These solutions 
are dual to interface conformal field theories at finite temperature. An approximate solution 
is first constructed using perturbation theory around a planar BTZ black hole. Numerical and 
exact solutions valid for all sets of parameters are then found and compared. Using the exact 
solution the thermodynamics of the system is analyzed. The entropy associated with the Janus 
black hole is calculated and it is found that the entropy of the black Janus is the sum of the 
undeformed black hole entropy and the entanglement entropy associated with the defect. 



1 Introduction 



The Janus solutions pQ in various supergravity theories provide interesting realizations of in- 
terface conformal field theories within the AdS/CFT correspondence [21 O H]. The simplest 
examples are constructed using an AdSd slicing of AdSd+i and making a massless field depen- 
dent on the slicing coordinate. Generically the scalar approaches two different constant values 
at the boundary of the space. Since a massless field is dual to a marginal operator the holo- 
graphic interpretation is that the coupling constant associated with the marginal operator jumps 
across the interface. The nontrivial profile of the massless scalar breaks the full SO(d,2) con- 
formal symmetry but due to the AdS^ factor the solution has SO(d — 1,2) interface conformal 
symmetr^ 

An important feature of the AdS/CFT correspondence is the fact that a black hole in the 
bulk of the AdS space is dual to a CFT at finite temperature [2]. It is natural to ask what 
is the bulk description of an interface CFT at finite temperature, i.e. a Janus black hole. In 
general this is a complicated problem due to the fact that a (planar) black hole in AdSj+i h as 
translational invariance in the d — 1 spatial directions. In a Janus solution the nontrivial scalar 
profile will break this translation invariance. Consequently a generic ansatz for a Janus black 
hole involves dependence on two spatial coordinates and hence solving nonlinear coupled partial 
differential equations, where numerical methods offer the only approach to solve the problem. 

In this paper we will focus on the simplest case of the Janus black hole in three dimensions. 
The simplification is two fold, both the BTZ black hole solution as well as the Janus solution 
are simpler than in AdS^+i with d > 2. Secondly we find an exact analytical solution^] for the 
Janus black hole which can be used to study analytically the physics of the interface CFT as 
well as test the numerical solution. Both simplifications are not present in higher dimensions 
but we believe that, as is often the case in AdS/CFT, there are valuable lessons to be learned 
from the lower dimensional exactly soluble system. 

The structure of the paper is as follows: In section [2] we give a brief review of the Janus 
solution (at zero temperature) and the holographic calculation of the entanglement entropy 
of the interface CFT. In section [3] a small dilaton perturbation with a Janus profile around 
a planar BTZ black hole is constructed and it is shown that the boundary stress tensor is 
unchanged to leading order in the perturbation. The numerical solution for the Janus black 
hole is discussed in section [4} The main technical problem, that the domain of the coordinates 
is not fixed but dependent on the solution, is solved by using the dilaton itself as a coordinate. 
In section [5] the exact solution is constructed and used to test the accuracy of the numerical 
solution. In section [6] we calculate various thermodynamic quantities associated with the black 
Janus solution. In particular we calculate the entropy, being the horizon length of the black Janus 
solution, using two methods. First, we use a boundary-horizon map based on bulk conservation 

1 Many generalizations of Janus solutions have been found over the years, see e.g. [5l l6l 171 151 l9l 1101 11111121 [13] . 
2 See also (131 H"5 16 for the other kinds of closely related exact solutions. 
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law considerations and second, we use a map which maps the boundary to the horizon using 
null geodesies. Interestingly both calculations give the same result and the entropy is the sum 
of the BTZ black hole entropy and the entanglement entropy of the zero temperature interface 
CFT. We close with a discussion of our results in section Various details of calculations are 
relegated to appendices. 



2 Janus system 

In this section we will briefly review the construction of the three dimensional Janus solution 
[H 113] and its holographic description as an interface conformal field theory (ICFT). 

The three dimensional Janus solution can be embedded into ten dimensional type IIB su- 
pergravity by the following ansatz for the metric 

ds 2 = e* /2 {dsl + ds 2 S3 ) + e~^ 2 ds 2 Mi (2.1) 

where M4 is either T 4 or K3. For the non supersymmetric Janus solution the dilaton and metric 
are independent of the coordinates of the three sphere and M4. All other type IIB supergravity 
fields are set to zero except the three- form field strength F3 |13| . whose explicit form is not 
needed in our later discussions. Upon dimensional reduction of the ten dimensional action on 
M4 x S3 one arrives at three dimensional action for the metric and <fi given b^] 

S = [ d^xJg (R - g ah d a 4>d h 4> + 2) (2.2) 



From this action the Einstein equation becomes 

Rab + 2g a b = d a 4>d b (t) (2.3) 
and the scalar equation of motion is given by 

da(V99 ab db^) = (2.4) 

For the static Janus solution the three dimensional metric and the dilaton take the following 
form 



ds 2 = dr 2 + f(r) 4> = 4>{r) (2.5) 



-dt 2 _ + d? 

e 

where we have set the AdS% radius to one for convenience. In [13] it was shown that the solution 



of the equations of motion (2.3) and (2.4) for this ansatz is given by 



/(r) = -(l + yi-2 7 2 cosh(2r)J (2.6) 

3 It is pointed out by Eoin O Colgain that this three dimensional Janus system can also be consistently embedded 
into the eleven dimensional supergravity using the extra dimensional geometry S2 x CY3 [17] . This is dual to an 
interface version of the chiral (4, 0) SCFT in two dimensions. 
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and 

W V2 V 1 + VI - 27 2 - v^T tanh(r) J K ' 

The Janus solution holographically realizes an ICFT where two CFTs defined on 1 + 1 dimen- 
sional half spaces are glued together over a + 1 dimensional interface. This can be seen as 



follows: The conformal boundary of the metric (2.5) has three components as r — > ±oo and 
finite £ > we can strip off the 1 /£ 2 factor and the boundary geometry are two copies of R x R + 
spanned by t, £. Note that the dilaton approaches two constant values at the boundaries 

hm <fr(r) = 00 + 4 lo S ( ^/PR^T^ ) = ^ ± M 

Without loss of generality we may set (fro = 0. For later use it is also convenient to express 7 in 
terms (fr as : 

7 = -j= tanh V2cfr as (2.9) 

The values of the dilaton on the boundary is dual to a modulus of the two dimensional CFT. 
The third boundary component is at £ = (i.e. the boundary of the AdS2 factor) This defines 
the interface where the two half planes are glued together. Hence the dual CFT is an ICFT 
where two CFTs defined on a half line are at different points in their moduli space. 

2.1 Entanglement entropy 

A useful observable in the ICFT is the entanglement entropy which is defined as follows. The 
space on which the CFT is living is divided into two domains A and £>. The total space of states 
H is expressed product H = ® Hg, where Hj\.,b is supported on A and B respectively. A 
reduced density matrix can be defined by tracing over all states in B, 

PA = tr He p (2.10) 

where p is the density matrix of the total system (at zero temperature this is just the projector 
on the ground state). The entanglement entropy associated with the domain A is then defined 
as 

Sa = -tTH A PAl°gPA ( 2 - n ) 

A holographic prescription to calculate the entanglement entropy in spaces which are asymp- 
totic to AdSd+i was presented in [181 [19]. F° r t ne domain A we denote the boundary dA which 
separates it from B. A static minimal surface which extends into the AdS^+i bulk and ends 
on dA as one approaches the boundary of AdS^+x- The holographic entanglement entropy can 
then be calculated as 

where Area(r_4) denotes the area of the minimal surface T_4 and G^ d+l ^ is the Newton constant 
for AdSd+i gravity. 
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In the present paper we consider the Janus deformation of AdS^ and the area A is an 
interval and the boundary dA are the two end points of the interval. The minimal surface is a 
space- like geodesic connecting the points. The geodesic which was used in [20] to compute the 
entanglement entropy chooses the £ coordinate as constant £ = £o> while r varies from — oo to 
+00. This corresponds to a symmetric region of width 2£o around the interface. 

The geodesic length is divergent and has to be regularized by introducing a cutoff e near the 
boundary [20j 

Area(r) = R AdSa J dr = R A dS 3 ( r oo(T) - r^T)) (2.13) 



where the regularized length can be read off from (2.6) 

r±oo = =F ( log e + - log - 2 7 2 - log(2£ )) (2.14) 

Hence 

Area(r) = Ra<IS 3 (^oo(r) - r_oo(r)) 

= iU dSa (2 log ^-log^l^)) (2.15) 

The holographic result has the same general form as the entanglement entropy calculated on 
the CFT side using the replica trick 



SU = |log- + log<M (2.16) 
6 e 

where we identify the length of the interval L = 2£o and e is the UV cutoff. The last term is 
boundary entropy (sometimes called g-factor [52]) which is associated with the degrees of freedom 
localized on the interface. In Rcf. [20j it was shown that in an expansion in the small deformation 
parameter 7 the holographic result ( |2.15 ) agrees to leading order with the weakly-coupled CFT 



calculation of the boundary entropy. Note that in the supersymmetric generalization of the 
Janus solution \'2'6\ |2"4] one finds exact agreement to all orders. In Appendix A, we redo the 
leading order computation using conformal perturbation theory, which is valid for the strongly 
coupled limit. 



3 Black Janus as a perturbation of the planar BTZ black hole 
3.1 Planar BTZ Black Holes 

We shall begin our discussion of three dimensional Janus black holes by studying the leading 
order corrections to the geometry and scalar field starting from the planar BTZ black hole 
solution. The planar BTZ black hole in three dimensions |25] can be written as 



ds z 



1 



(1 - z 2 )dr 2 + dx 2 + 



dz 2 



1 



(3.1) 
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where we take the range of coordinate x as (—00, 00). Of course the x direction may be com- 
pactified on a circle but we shall be concerned here only with the non compact case. The horizon 
is located at z = 1 and a convenient change of variables for checking its regularity is 



~2 -1 2 
Z = 1 — Z 



In this coordinate system the black hole looks like 

1 



ds' 



z 2 dr 2 + dx 2 + 



dz 2 



l-z 2 



l-z 2 

and now the AdS boundary is located at z = 1 while the horizon at z = 0. Since 

dz 2 



(3.2) 



(3.3) 



l-z 1 



+ z 2 dr 2 ~ dz 2 + z 2 dr 2 



(3.4) 



in the near horizon regime, there is no conical singularity if the euclidean time coordinate r is 
periodic with period 2ir. Therefore the corresponding temperature can be identified as 

1 



2vr 

The BTZ black hole with general temperature is described by the metric 

/2 



ds z 



Z' 



,-. 2 /2s , /2 . , ,2 . dz' 

(1 — a z Jar + dx + 



1 — a 2 z 



2 r /2 



which can be obtained by the scale coordinate transformation 



z = az t = ar x = ax 



from (3.1). The temperature for this scaled version of the black hole now becomes 

a 



T 



2tt 



(3.5) 



(3.6) 



(3.7) 



(3.8) 



Below we work mostly with the temperature T = (2ir) 1 and, using the above freedom of scale 
transformation, we shall recover the general temperature dependence whenever it is necessary. 



3.2 Linearized Black Janus 



Introducing a new coordinate y given by z = sin y, the planar black hole metric (3.1) can be 
rewritten as 



ds 1 



1 



sin 2 y 



cos 2 y dr 2 + dx 2 + dy" 



(3.9) 



Motivated by the form of this metric, we shall make the following ansatz for the black Janus 

^ 2 = ^^f + #^ <t> = <P(x,y) (3.10) 

A(x,y) B{x,y) 
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It is then straightforward to show that the equations of motion (2.3) and (2.4) reduce to 

0A) 2 - A d 2 A = 2A- A 2 0<j)) 2 (3.11) 
30B) 2 - 2B Pb = 8B 2 /A (3.12) 
dB- d4> -2B ^0 = (3.13) 



where we introduced the notation d = (d x , d y ). To the leading order, the scalar equation can be 
integrated as 

sinhx , o , , 

4> = 1 1 ^= + 0( 7 3 ) (3.14) 

J sinh 2 x + sin 2 y 

where we have used the Janus boundary condition 4>(x, 0) = 7 e(x)+0(7 3 ) with the sign function 
e(x). The leading perturbation of the metric part is of order 7 2 . For which we set 



with 



2 2 
A = Ao(l + ?-a(x, y) + 0( 7 4 )) B = B (l + \b{x, y) + 0( 7 4 ) 



Aq = sin 2 y Bq = tan 2 y 



The leading order equations for the metric part becomes 



2a — sin 2 y (Pa 



sin 4 y 



(sinh 2 x + sin 2 y) 2 



2 tan y d y b — sin 2 y Pb + 4a = 



(3.15) 
(3.16) 

(3.17) 
(3.18) 



With the Janus boundary condition whose detailed structure will be discussed later on, the 
solution can be found as a(x, y) = b(x, y) = q(x, y) where 



q{x,y) 



( sinh x 
V siny 



tan" 



_Y ( sinhx 
V siny 



+ 



sinh 2 x 
sinh 2 x + sin 2 y 



+ 2 + (ci - c 2 



sinhx 

I— 

smy 



(3.19) 



with ci and C2 being O(l) integration constanta^] Indeed checking that (3.19) solves eqs. (3.17 



3.18) is straightforward. Then the metric for the black Janus can be written as 

cos 2 y dr 2 + dx 2 + dy 2 ] + 0(~f A ) 



ds 



2 _ 1 - \i{x,y) 



sin 2 y 



Next we introduce a new angular coordinate fi that is defined by 

/ 7 2 \ sinh 2 

tan n + — ci = — 

V 4 / suit/ 



(3.20) 



(3.21) 



The above metric for the linearized black Janus can be written using the scale function /(/x) of 
the original Janus solution: Namely, the metric can be expressed in the following form 



ds 2 



f V 



-C-2 



sinh x + sin y 



cos 2 y dr 2 + dx 2 + dy 2 + 0{-i 4 



(3.22) 



4 Only the combination Ci — C2 is a true integration constant. We break it up to two since their roles are 
different in the discussions below. 
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where 



fifl) - sn 2 (K + (M + Mo)^ 2 ) ^ 



with 

4^(l±^/l- 2 7 2 ) (3.24) 
A; 2 = k 2 /4 = y + 0( 7 4 ) (3.25) 
Mo = i^(fe 2 )/^ + = | (l + g7 2 + 0(7 4 )) (3-26) 
To show this, we have used the expansion of the scale function /(//) given in the form 



As will be explicitly verified later on, the remaining part of the metric except the scale factor 
/ possesses a translational isometric direction along fx. Hence one may set C2 = without loss 
of generality. The zeroes of the function A and B occur at \i = ±/xo, which correspond to the 
boundary of the asymptotically AdS space. As a consequence the coordinate \i is ranging over 
[— /xoj Mo] as the case °f the original Janus solution. 

The choice of the integration constant c\ is simply related to the the choice of the coordinate 
patch of (t,x,y), which covers only part of the entire black hole geometry. Let us first consider 
the choice C\ = 0. Then the boundary occurs at 



smy 
sinhxl 



= tan (| - fio) + 0( 7 4 ) (3.28) 



This is then solved by 



37T 

siny(x) = -— 7 2 | sinhx| + 0(j 4 ) (3.29) 
16 



Therefore the validity of the coordinate is limited by 

3vr7 2 |sinha;|/16 < 1 (3.30) 

and, hence, the coordinate along the boundary become singular if 

\x\ > x cut ~ ln(2/ 7 2 ) (3.31) 

Only the region — f < /i < f can be free of any such coordinate problem. But this is simply 
a coordinate singularity, which may be removed by choosing a different coordinate chart. For 
instance let us consider the choice \c\ = | — hq + 0(7 4 ). For this case, one boundary can 
be solved by = for x > 0, whose solution is y = 0. For this side we do not have any 

coordinate problem but the other side of the boundary has again the coordinate singularity. The 
region specified by [Xq — ir < fi < fio of this coordinate chart does not involve any coordinate 
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problem. By the choice of \c\ = — § + /Uo + 0(7 4 ), the — /Uo < \i < — /Uo + tt region can be 
safely covered, which includes the other side of the boundary. 

For our further analysis of geometry below, we shall simply choose c\ = (together with 
C2 = 0), since the presence of the coordinate singularity can be ignored in the small 7 limit. 

3.3 Boundary stress tensor 

In this subsection, we shall construct the Fefferman-Graham metric to determine the boundary 
stress energy tensor. In order to use the prescription developed in Ref. [26], we introduce the 
metric in the following Fefferman-Graham form, 



d,Y 2 1 

ds 2 = -\ + -29^{X, X )dX^dX v (3.32) 



where (/i = 0, 1) denote the boundary coordinates^] and x = corresponds to the location 
of the boundary. In general one may expand g^ u by 

9flv (X, x) = 9$(X) + X 2 9$(X) + • • • (3.33) 

where gffi is the metric for the boundary system. In three dimensions, the boundary stress 
energy tensor is then given by [26] 

W = ^ " $(X) £l(X)g^(X)] + r, v {X) (3.34) 

where Tuu(X) is the scalar contribution for the stress energy tensor given by 

r^(X) = — [d^Bdvfo - y -f~ g^dafodpfo] (3.35) 

with <f>B denoting the boundary value of the scalar field. For our case, the boundary metric is 
given by 

g$ = diag(-l 5 1) = r? Mi , (3.36) 

since the boundary system is defined in the flat Minkowski space in two dimensions and the 
scalar contribution to the stress energy tensor vanishes since the scalar field is constant except 



X\ = 0. Let us first bring the metric in (3.22) to the form 



ds 2 = djq ( l _± c \_ dX 2 cot 2 Y ( l _± D \ (3 37) 

sin 2 Y sin 2 yV 4 > V 4 ) y ' 

where Xq = it. Introducing X(x,y) and y(x,y) by 

2 2 

X 1 (x,y)=x- 1 -X(x,y) + 0( 1 i ) Y(x,y) = y - ^y(x,y) + 0( 7 4 ) (3.38) 



5 We use this notation of boundary coordinates for this subsection only. In the subsequent (sub)sections we 
shall simply use t and x for the boundary coordinates for the notational simplicity. 
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and comparing the two forms of the metric to the leading order of 7 2 , one finds the differential 
equations, 



d y y -ycoty = q(x, y) d x y + d y X = 



(3.39) 



with 



C(x, y) = q(x, y) - d x X + y cot y 

y 

D(x,y) = q(x,y) + 

sin y cos y 



(3.40) 



The boundary conditions C(x, 0) = D(x, 0) = are required to have the standard form of the 
boundary metric rj^ v . The solution satisfying the boundary conditions is uniquely found by 

3 sinh 2 x + 2 



X(x,y) 



3 sin y cosh x tan 



-i 



sinh x 
siny 



+ (1 - cosy) 



sinh x cosh x 



cos y sinh x 



3 cosh 2 x + 1 
cosh 2 x 



tanh 



1 ^ - tanh" 1 A 



coshx 



coshx 



y(x,y) 



-3 cos y sinh x tan 1 — — — 3 sin y 



smy 



— siny- 



3sinh 2 x + 2r cosy i 
tanh tanh 



coshx 



1 



coshx 



coshx 



(3.41) 



with 



C(x,y) 



D(x,y) 



+ 



2 cosy 
cosh 3 x 

2(1 — cosy) 
sinh 2 x cosh 2 x 



, i cosy _i 1 

tanh tanh 

cosh x cosh x 



— secy- 



3 sinh 2 x + 2 



coshx 



_! cosy 
coshx 



tanh 1 — — tanh 1 



1 

coshx 



sin 2 y 



——2 . 2 + 3(1 - secy) 

sinh x + sin y 



(3.42) 



One may check that the location of boundary y(x) in ( |3.29 ) corresponds to Y = as expected. 

For the stress energy tensor, we now note that C(x,y) = D(x,y) = 0(y 4 ), which implies 
that the 0(j 2 ) terms have no contribution to the stress energy tensor. The rest is then straight- 
forward: Noting 



X = 2 tan 



Y 



one 



finds 



Therefore, one has 



y( 2 ) = ^diag(+l,+l) + 0( 7 4 ) 
' diag(+l,+l) + 0( 7 4 ; 



1 [jLV 



16vrG 



(3.43) 
(3.44) 
(3.45) 
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Hence to the leading order in 7, the stress energy tensor is independent of the deformation. 
Later on we shall show that T^ u is in fact 7 independent and the zeroth order result is all order 
exact. Finally recovering the temperature dependence by the scaling transformation, we have 

T^ = ^diag(+1,+1) (3.46) 
which agrees with that for the usual BTZ black hole. 



4 Black Janus at arbitrary 7 — Numerical Approach 

In order to obtain the form of the Janus black hole in the fully nonlinear regime of arbitrary 7 
we develop a numerical approach to solving the relevant Einstein+dilaton system of equations. 
It turns out that a successful implementation is surprisingly subtle, due to the special features 
of the Janus system. 

The most naive guess for the numerical ansatz would be 

' f ,y^V)+m^y) cog 2 y dT 2 + dx 2 + e V(x,y)-W( X , y ) dy 2l (4J 



ds 2 



sin 2 y 



with the dilaton given by (f)(x, y). This form of the metric ansatz makes it very easy to implement 
both constant temperature (which corresponds to Dirichlet boundary conditions for W(x, y) at 
y = tt/2 and Neumann for V(x, y)) and the Janus boundary condition. For numerics, we should 
map the infinite spatial coordinate into a finite interval e.g. by the mapping s = tanh x. However, 
this metric leads to several problems, some purely numerical and some, what is more dangerous, 
conceptual. Firstly, the discontinuous boundary condition for the dilaton with the jump at x = 
is very difficult to handle numerically. Secondly, it is far from clear what would be the domain 
of definition of the exact solution. The range of the bulk y coordinate is y € [0,7r/2], however 
the range of s is unknown. At the boundary the range is s € [—1,1], but in the bulk it may 
well be s S [— s max (y) ,s max (y)] with an a-priori unknown profile s max (y). This indeterminacy 
a-priori precludes any numerical treatment. 

In order to overcome the above difficulty, it is convenient to link the spatial coordinate to 
the value of the dilaton, since then the asymptotic range of the spatial coordinate is fixed by 
the definition of the Janus system, since at spatial infinities the dilaton is constant in the bulk 



and attains its asymptotic value ±4> as - This leads to the following ansatz (recall (2.8)) 



dr 2 

ds 2 = 2~ +A(y,s)dy 2 + 2B(y,s)dyds + C(y,s)ds 2 (4.2) 

tcin y 

</> = 4> as s (4.3) 

We have thus to deal with a non diagonal metric. Moreover it is not obvious what should be the 
boundary conditions characteristic of asymptotically AdS spacetime. To this end let us consider 
as a first approximation the BTZ black hole and use the linearized dilaton perturbation to fix 



the spatial coordinate according to (4.3). 
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The standard BTZ black hole metric is 



ds 2 



dr 2 



+ 



1 



tarn y sin y 



dy2 + JT^ 



(4.4) 



where we introduced s = tanhrc. Now the linearized Janus perturbation (3.14) takes the form 



s 2 + (1 — s 2 ) sin 2 y 

which leads to the following change of coordinates 

s sin y 



(4.5) 



\J\ — s 2 cos 2 y 

In these 'dilaton-adjusted' s — y coordinates the BTZ black hole takes the following form 



(4.6) 



ds' 



dr 2 



tan 2 y 



+ 



1 



1 — s 2 cos 2 y 



dy 2 2 s cos y ds dy 



ds 2 



sur 



y siny(l — s 2 ) (1 — s 2 ) 2 



which suggests the following ansatz for numerical computations 



ds 2 



dr 2 



tan 2 y 



+ 



1 



1 — s 2 cos 2 y 



K(y, s)dy 2 2L(y,s)ds dy M(y, s)ds 2 



sin 2 y 



siny(l 



1 



,2\2 



(4.7) 



(4.8) 



with very smooth solutions for (j) as = 0: K = M = 1 and L = scosy. Performing numerics with 
such an ansatz shows that even for very small 7 (equivalently (pas), the coefficient functions have 
different limits as y — > with s = 1 fixed and as s — > 1 keeping y = fixed. This leads to severe 
numerical problems and indicates that the pre-factor 

(4.9) 



1 — s 2 cos 2 y 



should be replaced by a suitable 7-dependent function. 

Fortunately, the Einstein-dilaton equations for the coefficient functions at y = reduce to 
ordinary differential equations with no y-derivatives, which can be solved exactly. Particularly 
relevant is the solution for K(0, s): 

f>as0-- s 2 )sinh v / 20as 



K(0,s) 



o 



with 



\/2(cosh V2~4>as — cosh y/2s) 
2 _ tanh \/2(j) as 

V2(j)as 

Taking into account the asymptotic BTZ metric at spatial infinity 



a 



K(y, 1) = M(y, 1) = 1 
we are led to modify the pre-factor 

1 

1 — s 2 cos 2 y 



L(y, 1) = cosy 



(4.10) 



(4.11) 



(4.12) 



(4.13) 
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to take into account these properties, namely we require that 

f(0,s) = (4.14) 
1 — s z 

sin y 

In addition it is convenient to have f(n/2,s) = 1 and dyf(n/2,s) = so as not to modify the 
form of boundary conditions at the horizon. A function which satisfies all the above properties 
is 



a 2 + (1 - q 2 )(1 - s 2 )sin 2 y 
1 — s 2 (l — a 2 sin 2 y) 

This leads us to the final ansatz for the numerical solution: 



dr 2 
tan 2 y 



ds2 = + f(v> s ) 



e K( y ,s) dy 2 2L(y,s)dsdy e M ^ds 2 
sin 2 y siny(l — s 2 ) (1 — s 2 ) 2 



(4.17) 



At the horizon y = tt/2 we impose the following boundary conditions: 

K(~,s) = H~,s) = d y M(~,s) = (4.18) 

The first condition ensures that the temperature is constant. At s = we impose boundary 
conditions following from symmetry 

d s K(y,0) = L(y,0) = d s M(y,0) = (4.19) 

4.1 Numerical details 



A feature of the Einstein's equations for the coordinate ansatz (4.17) is that we can pick three 
independent first order equations. These are the equation for <f>, the rr component of Einstein's 
equations E TT , and an appropriate linear combination of E yy and E ss . Subsequently we solve 
the above three equations numerically. 

We use a Chebyshev grid of = 20 or = 30 points in each dimension and use the PETSc 
library for solving nonlinear equations. We use a Python interface to the library (petsc4py). 
Unfortunately the convergence is poor so we had to use an automatic differentiation package 
(ADOL-C with Python bindings pyadolc) to compute the Jacobian and use LU linear solver 
from PETSc instead of the standard iterative ones. In addition the numerical solutions were 
found by gradually increasing the asymptotic value of the dilaton <p as by 0.05 from the BTZ 
value of (p as = and using the output from the previous value of 4> as as initial conditions for 
<p as + 0.05. In this way we generated the metric profiles up to (j) as = 10.0. 

The above numerical setup should be readily generalizable to the Janus black holes in higher 
number of dimensions for which, almost certainly an analytical solution does not exist. 
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Figure 1: On the left hand side we show e K ^ y ' s ^ for the numerical solution with 
on the right its relative deviation from the exact analytical solution. 



10.0, and 



5 Black Janus at arbitrary 7 — Exact solution in 3D 

As an alternative to the direct numerical solutions, one may attempt to use the above properties 



of the metric, such as the form of (4.10 ) to try to obtain an exact analytical solution. Remarkably 



this can be done for the three dimensional black Janus system considered in this paper. 

The exact Black Janus solution for arbitrary 7 turns out to be given by the following ana- 
lytical expression: 

du 2 



ds z 



cot udr + F(u, ip) 



sin 2 u 



4>l s f(<P)( 2 ■ 2 , 2 
I 7 sm u + cos u 



+ cot u (log f((p)) f dud(p 
cosh \/2(j) as ip sinh 2 \/2(f> as <p 



- + 



cosh \/2(f> as 2 cosh 2 \f2<pa 



dip' 



where 



and 



F(u,<p) 



f(<P) 



sm u + 



7 



cos 2 u 



cosh \[24>as 



(5-1) 



(5.2) 



(5.3) 



In this form we see the way that the exact metric incorporates the behavior (4.10). The similarity 



with the numerical ansatz made it easy to compare the above exact expression with the numerical 
solution for quite large value of 4> as = 10.0 with excellent agreement. In figure [T] we show the 



numerically obtained e K ( y ' s ^ coefficient function from (4.17) together with the relative deviation 



from the exact solution obtained from (5.1). 



It turns out, however, that passing to another coordinate system allows us to drastically 
simplify the metric. Indeed, let us change the u coordinate into w through the expression 



cot 2 « = f(ip)(w 2 - 1) 



(5.4) 
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Then the exact metric transforms into 



ds 2 = f(<p) 



[w 2 - l)dr 2 + 



duu 



+ 



6 2 
as 



w 2 — 1 7 



(5.5) 



This can be further simplified introducing the coordinate p instead of w through w = cosh p. 
With this substitution, we obtain a remarkably simple form of the exact Janus black hole metric: 



ds z 



f(<p) sinh 2 pdr 2 + dp 2 + ^-f(<p)d<p 2 



V 



/(/i) sinh 2 pcir 2 + dp 2 + d\jf 



where we have used the fact (jj^j V~f(p) = ll4>as with /(//) given by (3.23). 



(5.6) 



6 Thermodynamics and Entropy 
6.1 Thermodynamic quantities 

In this section, we shall be describing thermodynamic properties of our Janus system based on 
the solution of the previous section. First we turn to the stress energy tensor. We have already 
discussed the stress energy tensor for the case of the linearized black Janus and claimed that 



the expression in (3.46) is in fact exact to all orders in 7. Let us argue this point first. From 



the exact metric we note that the deformation parameter 7 dependence occurs only through the 
function /(/i). In the near boundary region of fj, ~ ±/Uo, the scale function /(/x) can be expanded 
as 



1 



1 



!+o(m-Mo) +&l(7)(m-Mo) + MtXm - fJ-o) 



(6.1) 



and 7 dependence can only appear in the higher order coefficients &2n with n > 2. From 
dimensional analysis it follows that (/jl — Ho) 2n ~ \ 2n /x 2n in the near boundary region where x 



is the Feffer man- Graham coordinate introduced in (3.32). Hence gj$ and gffi are 7 independent 



since its dependence enters only at higher orders of x- Therefore we conclude that the stress 



energy tensor in (3.46) is all order exact. 



One can reach the same conclusions purely from the field theory perspective. Consider the 
energy- momentum conservation equations: 



d t T a + d x T xt = 



d t T xt + d x T x 







(6.2) 



Due to time reversal invariance we have T xt = 0. Then T xx is x independent and, due to 
tracelessness, so is T u = 0. So T^ v is diagonal and constant in space and time. For an infinite 
system that we consider explicitly here, the value of the energy density equals the value of the 
energy density in the asymptotic region and so equals the BTZ value. 

Of course our expression of stress tensor is strictly valid only when the size of the boundary 
system L goes to infinity. Dealing with the finite size system, which in general involves the finite 
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size effect, is not a simple matter. According to the recent proposal [27J, the boundary dB of 
a boundary system, can be dealt with in a rather simple manner by introducing corresponding 
bulk boundaries dM. This hypersurface is extended into the bulk in an appropriate manner 
from dB and introduces an extra contribution to the thermodynamic quantities. This formalism, 
however, is developed only for boundary conformal field theories and it is not clear whether it 
is applicable to our case or not. Thus in this note we shall be only concerning about the limit 
where the system size L goes to infinity. 

In the next two subsections, we shall be calculating the entropy of the black Janus system 
utilizing two different methods. Since the translational symmetry in the x direction is broken 
by the Janus deformation, the entropy density should be position dependent. The entropy is 
a quantity defined at the horizon, whereas the interface CFT and its stress tensor is naturally 
defined at the boundary side of the geometry. Hence in order to talk about the entropy, one has 
to relate the horizon side to the boundary side. Especially one needs a map which relates the 
boundary coordinate x to the horizon coordinate in order to use the data found at the horizon 
of the black hole. Below we shall discuss two methods for the boundary-horizon map, which will 
be used to determine the entropy of the boundary system. We shall find that the both methods 
lead to the desired entropy in the large L limit to all orders in 7. Using the map, one may wonder 
whether one can define an entropy density by evaluating the corresponding horizon length scale 
divided by AG. The resulting expressions for the entropy density by the two methods turn out 
not to agree with each other even in the large L limit. This, however, is not a problem. As 
discussed in [28] , the entropy density alone is not fully well defined in a gauge invariant manner 
even in the semiclassical limit. We shall get back to this issue later in the discussion section. 

6.2 Method 1 

As stated previously, the entropy is a quantity which is defined at the horizon whereas the 
interface conformal field theory is defined at the boundary of the geometry. The stress energy 
tensor is also determined by the behavior of the geometry near the boundary. The question is 
how to connect a boundary region given by AB to the corresponding horizon region A\B\ as 
depicted in Figure [2] In other words the question is how to find the two bulk lines AA\ and 
BB\ which join the image at the horizon to the boundary region (of a coordinate size L). This 
may be done as follows: Consider the conserved current |29j 



where £ c is the time translation Killing vector satisfying the Killing equation V £& + Vf,£ a = 0. 
In our case, £ a = 5 a0 and one may find the boundary lines from the condition 



with t = constant. This basically says that there is no contribution from the bulk lines for the 
integral of Q a over the region surrounded by AAiB±B. This condition then ensures the first 




(6.3) 



dx a Q a = 



(6.4) 
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law of thermodynamics that connects the change of entropy to the energy and the length of the 
system. 




Figure 2: The finite temperature Janus system, whose spatial extent is given by the segment AB, is 
corresponding to the bulk region AA\B\A. The dotted lines describe the identification for the undeformed 
geometry with 7 = 0. 



For the analysis of the exact geometry, we shall use the metric in (5.6). As in the previous 



section, we shall consider the boundary system of spatial size L with the interface at the cen- 
ter. (Namely x £ [-L/2, L/2] with x as the standard boundary coordinate.) Following above 
prescription, the corresponding bulk lines satisfy 



d<p _ g pp d 9 g tt 
dp g w dpgu ' 



This is solved by 



tanh 



tanh 



) = ( COSll P ^ COBh^ a 

(<t>aa\ \ 



cosh pa 



) 



For the line starting from xa = L/2, 



1 , L 

= tanh — 

cosh pa 2 



(6.5) 



(6.6) 



(6.7) 



This identification of the boundary coordinate for exact black Janus solution is carried out in 
Appendix C. The corresponding horizon point A\ is described by 



tanh ( < t >as< ^ Al \ = tanh f tanh — 

V y/2 J V A /9/ V 2 



V2' 



V2<t> a 



(6.8) 



together with pa x = 0. This is exactly the boundary-horizon map obtained by the present 
method. Then the horizon length Lh can be obtained by 



Lh 
2 



where we used 



pa-, 1 
f((p)dtp = - log 
7 Jo 2 



cosh V2(f) a 



l + (tanhfj V 
l-(tanhfj V 



l-2 7 2 



(6.9) 



3.10) 
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The entropy then becomes 



SL = AG l ° g 



1 + ( tanh § 



V1-27 2 



1 - ( tanh \ 



By expanding in 7 2 , one may find the leading order behavior as 



1~ sinh L log tanh 2 \ + 0( 7 4 ) 

oGr 2 



For the large size limit the finite part of (|6.11|) becomes 

AS L 



1 



00- 4G log(l/ V /l- 2 7 2 ) 



AG 



+ 0( 7 4 



3.11) 



.12) 



3.13) 



This expression precisely coincides with the result from the method of the holographic entan- 
glement entropy. 



6.3 Method 2: a boundary horizon map based on null geodesies 

In this subsection, we shall carry out the boundary horizon map based on the null geodesies 
emanating from the boundary in a hypersurface orthogonal manner. This construction is based 
on the concept of the light-sheet of the holography [30] and widely used in application of fluid- 
geometry correspondence [281 ED E21 [33] • 

Since null geodesies are the same for conformally equivalent metrics it is enough to determine 
them for the metric 



where we have 



The geodesic equations read 



and the null condition is 



ds 2 = - sinh 2 pdt 2 + dp 2 + dfi 2 



sinhp 



cot u 

77W) 



7 



f(s)ds 



t" + 2 coth p p't' 
p" + cosh p sinh p t' 2 



fx + p = sinh p t 









(6.14) 

(6.15) 
(6.16) 

(6.17) 
(6.18) 
(6.19) 

(6.20) 
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It is convenient to parameterize the geodesies with p instead of the affine parameter A. One can 
check that the solution is 

ci coshp 

t(p) = C2 + arctanh (6.21) 



c\ — sinh p 
coshp 

/j,(P) = arctan y + C3 (6.22) 

c 2 — sinh 2 p 

I c? — sinh 2 p 

dp/dX = ^ — (6.23) 

smhp 

As mentioned above, p at the horizon is 0, while approaching the boundary point leads to 

sinh 2 p A = — 3_ ( 6 .24) 
sinh x 



in agreement with (6.7). In order to determine which geodesic to take we have to fix the direction 
of the null geodesic emitted from the boundary. In the following we utilize the (u, to) coordinates 
which are defined in appendix [Cj Near the boundary this requirement means that we follow a 
curve of fixed uj. We have to translate this condition to a condition on the derivative 

f (6.25) 
dp 



at p = pa- To this end we have 



djs = d^ds = <K 1 rj^ds 
dp ds dp 7 v dp 



where the derivatives have to be taken along constant uj. We can evaluate the last derivative by 
differentiating 

S inhp= C ° t( ^ (6.27) 

v/0) 

with respect to p and taking the limit s — > 1. From this we see that ds/dp is finite which means 
that 

dfI ' 00 (6.28) 



dp 

Evaluating ^ gives 

dfj, s'mhp 



3.29) 



dp ycf — sinh 2 p 
Hence this condition enables us to identify the constant 

ci = sinh pa (6.30) 

Now the geodesic takes the following form in the jj,-p plane: 

COShp 7T 

fi = arctan = — — + ^0 (6.31) 

'sinh 2 pa — sinh 2 p 
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where /Uo corresponds to cj) = 1. Evaluating the above expression at the horizon we get 

1 7T 7T 

= arctan — h = arctan sinh x h fio 

smh 2 2 



3.32) 



This is the boundary- horizon map based on null geodesies. By inspection we see that it is 
different from the one in ( |6.8| ). 

If we take x — > oo (oo — > oo), and we recover \in — > ^$ which is the expected result. However 
one may be puzzled by the opposite limit. When oj — > which should correspond to x — > 0, [in 
does not approach but rather 



Mo 



7T 



(6.33) 



As an aside, from the general formula for dp/dfi we see that at the horizon the geodesic is always 
perpendicular to the horizon so the latter condition does not allow us to discriminate between 
geodesies. 



Based on the above boundary- horizon map in (|6.32), the horizon length is given by 

1 + sn(K +/ UH-, k 2 ) 



L H = 2 £ H dny/fM = log 



1 — sn(K + fiH, k 2 ) 



We also rewrite (6.32) as 



sin(^o — ^h) = 1/ coshxA = 1/ cosh 



L 



3.34) 



(6.35) 



From this formula, we would like to identify the leading order correction of the entropy. Noting 



7T 



MO 



1 + o7 2 +0(7 



one finds 



sin hh = tanh xa + 



3vr 7 2 



+ ... 



16 cosh xa 

In addition, the Jacobi sine function can be expanded as 

sn(K + fi H , k 2 ) = sin + — ^ - 3/i# cos + sm^ H cos 2 (ihj + 
From this, one finds 



L H = L + 



Interestingly, this leads to 



7 



L 



(A 

V msn 



tanh h 3 cosh — arcsin , 

2 2 V cos h L 



0( 7 4 



(6.36) 
(6.37) 

(6.38) 

(6.39) 
(6.40) 



as L becomes infinity. Thus we obtained the expected result to this order. Let us now rewrite 

dn(re + (/x - Mff), k 2 ) + cn(K+(^ - Hh), k 2 ) 



(6.34) as 



L H = log 



dn(K + (no - Vh), k 2 ) - cn(K+(^ - Vh), k 2 ) 



3.41 
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Then using 

Mo — = arcsin (seen— ^ (6.42) 

one finds 

L H ^L + log-J== (6.43) 
V 1 - 2 7 2 

Again this leads to 

AS i _ = i Llog 7r =L= (6.44) 
which agrees with the result of the holographic entanglement. 

6.4 Remarks on the agreement between the two methods 

Let us discuss now what aspects of the boundary-horizon map are probed by the above calcula- 
tions of the agreement between the entanglement entropy and the two computations using the 
two choices of the boundary- horizon map. 

If the size of the system would be finite (with a circle compactification) , then the total 
entropy would be obviously completely independent of the choice of the boundary-horizon map 
as it has a purely geometrical definition as the area (length in the 3D case) of the horizon. In the 
case of the infinite system that we consider, we have to subtract off the extensive contribution - 
and hence a potential difference may arise only in the upper limit of integration of the horizon 
area element. Therefore we are probing differences between the boundary- horizon maps in the 
asymptotic near BTZ region. Unfortunately, so far we do not have a test which would be 
sensitive to the finer details of the boundary-horizon map closer to the defect. 

6.5 First law of thermodynamics 

From the previous investigation, the entropy of the system in the large size limit is given by 

ttTT 

S = ^G +Si (6 - 45) 



where Sj denotes A5/ /=00 in (6.44) and we have recovered the temperature dependence. Sj is 



temperature independent and can be interpreted as the localized interface contribution to the 



entropy. From the stress tensor in (3.46), the energy and pressure can be identified as 



„ itT 2 L ttT 2 . A . 

E = ^r p =4o i6A6) 

Therefore one can check that the first law of thermodynamics 

TdS = dE+ pdL (6.47) 

holds precisely for our Janus system. Of course one should note that our investigation is valid 
only for the large size limit. 
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7 Conclusion 



In the present paper we have considered a supergravity dual to a three-dimensional interface 
conformal field theory at finite temperature. The supergravity fields which are turned on are the 
metric and the dilaton. The interface is realized as the boundary between two domains (x > 
and x < 0) with differing values of the vacuum expectation value of the operator dual to the 
dilaton ±(j) as . 

The undeformed finite temperature case corresponds to the well-known BTZ black hole. We 
have started from a linearized analysis of the scalar perturbation with small 4> as of the BTZ 
black hole building up intuitions concerning the general structure of the solution including the 
extraction of the boundary energy-momentum tensor which is quite intricate in coordinates 
natural for the Janus solution. 

Using this knowledge we have formulated a scheme suitable for the numerical computation of 
the exact Janus black hole for arbitrary <j) as , the key obstacle being an a-priori lack of knowledge 
about the coordinate domain of the exact nonlinear solution. We overcame this problem by 
linking a spatial coordinate to the value of the dilaton. 

In the case of three-dimensional system we found an exact analytical solution of the finite- 
temperature Janus black hole. Let us emphasize that, due to the coupling between gravity and 
the scalar field, the three-dimensional Einstein-dilaton system is nontrivial in contrast to pure 
three-dimensional gravity. The existence of an analytical exact solution was in fact completely 
unexpected for us. 

This exact solution is very interesting from various points of view. On the general relativity 
side it provides an example of a black hole in equilibrium with a nonuniform horizon. It would be 
very interesting to explore such features as Hawking radiation and temperature in this setting. 

From the AdS/CFT perspective, such a black- hole gives a dual description of a three- 
dimensional interface CFT at finite temperature, which may have condensed matter applications. 
On a more theoretical side, such a nonuniform black hole provides a theoretical laboratory for 
investigating various issues dealing with entropy of the dual field theory. In particular, we can 
probe various maps between boundary points and horizon points which have been proposed in 
the literature, and critically examine the problem whether local entropy density can be defined 
at all. 

As a step in this direction we have evaluated the total entropy deviation from the BTZ an- 
swer using two different boundary-bulk map prescriptions finding agreement with entanglement 
entropy calculations. 

There are numerous directions for further research. First, one may also consider ICFT 
defined on a circle with two interface points. Of course there is a corresponding Janus solution 
dual to this compact version of ICFT. This system can serve as an ideal setup to study finite size 
effects of a finite temperature ICFT. Secondly one can investigate the gravitational aspects of the 
nonuniform horizon, generalize the setup to higher number of dimensions (where probably one 
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would have to rely on numerical methods). One new nontrivial feature in higher dimensions may 
be the appearance of nonuniform energy density which is kinematically ruled out in the three- 
dimensional case. It would also be an interesting question to explore if and how the presence of 
an interface modifies the occurrence of a phase transition (such as a confinement /deconfinement 
transition) in AdS/CFT. 
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A Computation of entropy correction based on the conformal 
perturbation theory 

In this appendix we are interested in computing free energy defined by 

/3F = -logtre^ (A.l) 
Perturbing the above from Hq of BTZ system, one has 

(3F = f3F + 7 2 /3F 2 + 0( 7 4 ) (A.2) 
The 0(7) contribution vanishes because 

(C(-iT,x)) = (A.3) 

where the expectation value is evaluated with respect to the BTZ system with Hq. And one can 
find 

1 p/3 r/3 poo poo 

f3F 2 = - dr dr' dxe(x) dx'e(x')(£(-iT, x)L{-ir' , x')) (A.4) 

2 jo Jo J -00 J -00 

where e(x) is the sign function. Note that the two point function [15] is given by 

(C(-iT,x)£(-iT J ,xf)) = —-2 = 1 2 ( A - 5 ) 

lb7r G [cos(r — t') — cosh(x — x') — ie] 



for the BTZ background with (3 = 2ir. Using the integral 

Jo (cosx + Q) 2 " 27F (Q2 _ 1)3 /2 

/3F2 can be rearranged as 



2n 1 O 

dx 7 T^v? = 2?r (A.6) 



1 f°° , f°° , / / cosh(x - x') + q 2 , k . 

BF 2 = —r= dx dx'e(x') ^ r (A.7) 

iGJo J -00 V ' ((cosh(x - x') + q 2 ) 2 - I) 3 / 2 
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with q = it. Using the symmetry of the integrand, the integral may be further arranged as 

1 f L/2 dx f x d > coshx' + q 2 
P 2 ~ ~2GJo X Jo X ((coshx' + g 2 ) 2 - I) 3 / 2 ( ' ' 

where we also introduce a finite system size L. This regularized integral is finite and has the 
expansion, 

(3F 2 = c_i4 + c_i- + c (L) + 0(q) (A.9) 
q A 2 q 

where L dependence in the leading term is introduced by the dimensional ground. c_i is not 



physical since it reflects our choice of the vacuum energy level. We evaluate the integral (A.8) 
by the replacement of the integrand by 

1 fL/2 fx (1 + ,7 2 )coshx / 

£F 2 = / dx / dx' ^—^r 1 ^ A.IO 

2G7 Jo ((l+g 2 )sinhV + 2g 2 + g 4 ) 3/2 

which does not affect c_i, c_i and Co terms in the expansion. This fact may be checked 

2 

numerically, whose details are in Appendix B. Thus, by the replacement, the physical term 
cq(LT) is not affected. 

Then the integral may be carried out explicitly leading to 



1 l + g 2 (g2 + 1 ) cogh L + /( g 2 + 1 ) cosh 2L +g 4 + ^2_ 1 

^ = 2ovT7 log rr^TTgr? (A ' U) 



which can be expanded as 



2G 

Therefore one has 



AS=^+0( 7 4 ) (A.13) 



as the size L becomes large. 

Note that our computation here does not care about possible finite size effect. Namely any 
effect of boundary conditions at x = ±L/2 was not incorporated into the computation above. 
Including the finite size effects on the field theory side is an interesting open problem. 

B Equivalence of and T% to 0(q°) 

We are interested in the integral 

fL/2 rx 

J(q,L)= dx dx'l(q,x') (B.l) 



JO 



where 

-w \ coshx + o 2 , 

l(q,x) = ^ (B.2) 

((coshx + g 2 ) — l) 3 / 2 
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We will show that 



J(q,L) 



L 



+ * coth^ +0(q) 



V2q ' 2"~~ '2 

Using an integration by parts and a change of variables, one may rearrange J(q, L) as 



J(q,L) 



dx 



L 

2 



l(q,x) = q 



dy 



L 

2 



qy 



i(q,qy) 



(B.3) 



(B.4) 



We then expand the integrand into a power series in q and perform the integral order by order. 
And then expand the result again into a power series in q, which leads to the series expansion 
of (B.3) in L. For instance by expanding the integrand to the order q G , we get 

L 5 L 7 



J(q,L) 



L 



1 



1 L_ ^ i ^ _ ^ , 
4<? 2 ^2g + £ 12 720 + 30240 1209600 + ^ 



(B.5) 



which agrees with the L expansion of (B.3). This check can be pushed further to the higher 
orders in L. 



C The boundary metric and coordinates of exact solution 



In this appendix we shall construct the boundary metric from the exact solution (5.1). For this 
we have to pass to the Fefferman- Graham coordinates but just at the leading order. In order to 
approach the boundary we will take the limit u — > in a correlated manner with ip — >■ 1 namely 

<p = l--u 2 (C.l) 



keeping u fixed. From the exact metric, one finds 

2<f>as 

7 sinh 2 p 

The exact Janus BH metric reads then in the (it, oo) coordinates at leading order in u: 



w = . 2 (C.2) 



where 



ds^ = d 4 + d 4 + 11 ^ (C.3) 



H(u) = ^— (C.4) 



The boundary coordinate is then found from 



<H(uj)du} = dx (C.5) 
Straightforward integration leads to 

X = log (y^U + y/jUJ + 24> as ^j (C.6) 

which can be inverted to yield 

oj = as sinh 2 x (C7) 

7 

Finally using (C.2), one finds 

— \- = sinh 2 x (C.8) 
sinh p 
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